Van Hove Singularity and Apparent Anisotropy 
in the Electron-Phonon Interaction in Graphene 
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We show that the electron-phonon coupling strength obtained from the slopes of the electronic 
energy vs. wavevector dispersion relations, as often done in analyzing angle-resolved photoemission 
data, can differ substantially from the actual electron-phonon coupling strength due to the curvature 
of the bare electronic bands. This effect becomes particularly important when the Fermi level is close 
to a van Hove singularity. By performing ab initio calculations on doped graphene we demonstrate 
that, while the apparent strength obtained from the slopes of experimental photoemission data is 
highly anisotropic, the angular dependence of the actual electron-phonon coupling strength in this 
material is negligible. 



The energies and lifetimes of charge carriers in solids 
are significantly affected by many-body interactions in- 
cluding those with electron-hole pairs, plasmons, and 
phonons. Angle-resolved photoemission spectroscopy has 
emerged as an ideal tool for directly probing the effects of 
these interactions on the electron quasiparticle dynamics 
with good energy and momentum resolution 

In particular, the low-energy electron dynamics at 
metal surfaces [lin layered materials, such as mag- 
nesium diboride Q, graphite 0], and cuprate super- 
conductorsjlj ,and in single layer graphene 0, [g, 0, 
i, S M lUT El El HI M El, is significantly af- 
fected by the electron-phonon interaction. Since the 
electron-phonon interaction generally manifests itself as 
a kink in the quasiparticle dispersion relations measured 
by angle-resolved photoemission spectroscopy [13, EH, 
it is common practice to determine the strength of the 
electron-phonon coupling by taking the ratio between the 
group velocity at the Fermi level and below the phonon- 
induced kink [U, [n| [2(| [2l| . In the cases where this sim- 
ple procedure is not applicable, more complicated self- 
consistent algorithms 0, [H| become necessary to ex- 
tract the electron-phonon coupling strength. However, 
the application of these methods requires the knowledge 
of several adjustable parameters and is subject to some 
arbitrariness. Therefore, assessing in the first instance 
the validity of extraction procedures based on the linear 
slopes of the photoemission data is an important issue. 

Graphene @,[25|,[2i|[2Zl,[ii] is an ideal system to inves- 
tigate these effects. Indeed, the Fermi level of graphene 
can be tuned over a wide energy range by chemical dop- 
ing [1, [29[ or by gating [2(| [271 [3(| , and can almost be 
aligned with the van Hove singularity at the M point of 
the two-dimensional Brillouin zone [8[. 

In this work, we show that the apparent electron- 
phonon coupling strength in doped graphene obtained 



from the linear slopes of the renormalizcd quasiparticle 
dispersions, as calculated from first principles, is highly 
anisotropic, in good agreement with experimental re- 
sults @. On the other hand, the phonon-induced elec- 
tron self-energy is found to be only weakly dependent on 
the wavevector in the Brillouin zone. As a consequence, 
the actual electron-phonon coupling strength is isotropic. 
The apparent anisotropy of the electron-phonon interac- 
tion is shown to arise from the curvature of the bare 
electronic bands of graphene, which is strongly enhanced 
in proximity of the van Hove singularity at the M point. 
Our findings are relevant to the interpretation of pho- 
toemission spectra in materials where the Fermi level is 
aligned with a van Hove singularity, such as the hole- 
doped cuprates at optimal doping. 

The mass-enhancement parameter or electron-phonon 
coupling strength Ak of an electronic state with wavevec- 
tor k on the Fermi surface can be expressed through the 
energy derivative of the real part of the self-energy arising 
from the electron-phonon interaction [3l| : 
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E-p being the Fermi level. Within the Migdal approxima- 
tion, which corresponds to considering the non-crossing 
electron-phonon self-energy diagrams, this quantity can 
be calculated by [13, HI: 
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where e„k is the energy of an electron in the band 



with wavevector k, 



the energy of a phonon in the 
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branch v with wavevector q, and / n k and the 
Fermi-Dirac and Bose-Einstein occupations, respectively. 
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The integration is performed within the two-dimensional 
Brillouin zone of area A^z ■ The electron-phonon matrix 
element <? TOn) „(k, q) = (mk + q|AV4j„|nk) is the ampli- 
tude for the transition from an electronic state |nk) to 
another state |mk + q) induced by the change in the self- 
consistent potential AV^„ generated by the phonon \qu). 
The technical details of the calculations are reported in 
Ref . 0. 

It can be shown [l7| that the actual electron-phonon 
coupling strength in Eq. (fl]) can also be written as 



"£(£f) 

v k (E F ) 



1 



(3) 



where v k (Ep) and v k (Ep) are the renormalized and the 
bare velocities at the Fermi level, respectively. While 
Eq. ([3]) is useful for theoretical analyses, it cannot be 
used directly in determining the electron-phonon cou- 
pling strength from the experimental data since the bare 
group velocity is merely a conceptual tool and cannot 
be measured. To circumvent this difficulty, from the ex- 
perimentally measured low-energy photocmission spec- 
trum, the electron-phonon coupling strength is usually 
extracted [l|,l8| by taking the ratio of the renormalized ve- 
locity below and above the phonon kink. This procedure 
rests on the assumptions that (i) well beyond the phonon 
energy scale the bare velocity is fully recovered and (ii) 
the bare band is linear over the energy range considered. 
We denote the value obtained from this procedure as the 
apparent electron-phonon coupling strength: 
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where AE is taken slightly larger than the phonon energy 
Wph so that the energy E = E F — AE falls below the 
phonon kink. 

In Fig. [1] we compare the apparent coupling strength 
A^ pp obtained from our first-principles calculations and 
that extracted from the experimental photoemission 
spectra of graphene at four different levels of doping [1| . 
To determine v k (E F — AE) in Eq. (T5]) from our first- 
principles calculations, we considered the slope of the 
quasiparticle band at the energy AE =0.3 eV below the 
Fermi level (33[. We have checked that as the energy 
range AE varies within the interval 0.2~0.4 eV, the ap- 
parent strength A app changes by less than 10 % along 
both the KM and KT directions. Theory and experiment 
are in good agreement with each other for all doping lev- 
els considered. For graphene at the highest doping level 
(E F — Ep, = 1-55 eV, where Ep, being the energy at the 
Dirac point) the agreement between theory and experi- 
ment is slightly worse along the KT direction. 

Figure Ufa) shows the calculated apparent strength 
A app [Eq. (H}] as a function of doping for two different di- 
rections in the Brillouin zone of graphene. The apparent 
electron-phonon coupling strength is highly anisotropic 
and can become as large as A app = 2 along the KM di- 
rection for the doping levels considered here. 



Now we consider the actual electron-phonon coupling 
strength A as obtained from Eq. © [Fig. [5^b)]. The ac- 
tual strength increases monotonically with doping, reach- 
ing A = 0.22 when E F — Ep, = 1.5 eV. At variance 
with the apparent strength, the actual strength A does 
not depend on the direction of the wavevector k. We 
have checked that this holds for any path through the 
K point. The present results indicate that the actual 
electron-phonon coupling strength in doped graphene is 
extremely isotropic. Thus, the actual strength can differ 
substantially from the apparent strength, the more so as 
the Fermi surface approaches the van Hove singularity at 
the M point. 

In order to analyze in detail the angular dependence 
of the electron-phonon coupling in graphene, we cal- 
culated the phonon-induced electron self-energy using 
Eq. (1) of Ref. [H While in Ref. [U the electron self- 
energy was evaluated under the constraint E — e k , we 
here consider the complete energy-dependent self-energy 
J^k(E). Figure [3] (a) and (b) show the real and the imagi- 
nary part of the electron self-energy in n-doped graphene 
(Ep — Ep, = 0.64 eV). The wavevector is varied along two 
different paths indicated in the upper left corner of Fig. [3] 
The dependence of the self-energy on the wavevector k is 
found to be extremely weak, the variation along the path 
considered in Fig.[3]being less than 3 meV for a given en- 
ergy E. The insensitivity of the electron-phonon coupling 
strength Ak to the wavevector k [see Fig. [2jb)] is fully 
consistent with the finding on the self-energy. Figure [3J 
also shows that, while the self-energy is highly isotropic, 
the corresponding spectral function exhibits significant 
angular dependence due to the anisotropic dispersion of 
the energy bands in graphene. 

Recently, the observed anisotropy § in the appar- 
ent electron-phonon coupling strength has been related 
to foreign atoms based on calculations of CaC6 layers 
with the dopants arranged periodically in atop sites on 
the graphene plane [14} . Our calculations clearly show 
that the anisotropy in the apparent strength A app is al- 
ready present without invoking the possible effect of the 
dopants. 

Our investigation of doped graphene allows us to 
discuss some general aspects of the extraction of the 
electron-phonon coupling parameters from angle-resolved 
photoemission data. By expanding the energy depen- 
dence of the velocity to first order, we can rewrite approx- 
imately the apparent electron-phonon coupling strength 
in Eq. g]) as: 
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In Eq. $5§, v'^(Ep) is the energy derivative of the bare 
velocity, and we assumed that well below the phonon 
kink the bare and the renormalized velocities coincide. 
Equation ([5]) shows that, whenever the band velocity de- 
creases with decreasing binding energy (i.e., v k < 0) as 
is the case for graphene along the KM direction, the ap- 
parent electron-phonon coupling strength always exceeds 
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FIG. 1: (color online) Polar plots of the apparent electron-phonon coupling strength A^ pp on the Fermi surface around the K 
point in the Brillouin zone. Filled red squares and empty blue circles represent results from the ab initio calculation and from 
the experimental photoemission spectra, respectively. The lines are included as a guide to the eye. Different panels correspond 
to different doping levels. Note that the scale in (d) is different from that in (a)-(c). 
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FIG. 2: (color online) The apparent strength A app (a) and 
the actual electron-phonon coupling strength A (b) calculated 
along two different directions in the two-dimensional Brillouin 
zone of graphene: along KM (solid red lines) and along KT 
(dashed blue lines). Along the KF direction, A app can even 
become negative [cf. discussion around Eq. ([5}]. 
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FIG. 3: Calculated (a) real and (b) imaginary part of the elec- 
tron self-energy E(_E, k) and (c) logarithm of the correspond- 
ing spectral function A(E, k) arising from the electron-phonon 
interaction in n-doped graphene (Ef — Er> = 0.64 eV), along 
two different directions KM and KF in the Brillouin zone. 



the actual strength. 

To illustrate this point, we consider in Fig.|¥]the quasi- 
particle bandstructure for a model system obtained by as- 
suming an Einstein phonon spectrum with phonon energy 
w p h = 0.2 eV, a constant density of states near the Fermi 
level, and a constant electron-phonon coupling strength 
A=0.2. Within this model, the real part of the electron 
self-energy due to electron-phonon interaction reads [l8| : 



Re £(£) = 
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having included broadening F = 10 meV for convenience. 
As shown in Fig. BJa), the apparent strength A app = 0.21 
constitutes a good approximation to the actual strength 
A = 0.2 when the slope of the bare electronic band does 
not change appreciably within the phonon energy scale, 



i.e. v'£(Ef)AE <C v-k(E-p). However, in the case where 
the bare velocity decreases with decreasing binding en- 
ergy, the apparent strength A app = 1.8 differs signifi- 
cantly from the actual electron-phonon coupling strength 
A = 0.2, consistent with Eq. © [Fig. HJb)]. In the lim- 
iting situation where the Fermi level is aligned with the 
van Hove singularity (as in heavily doped graphene), the 
velocity at the Fermi level vanishes while the velocity be- 
low the phonon kink energy remains finite. As a result, 
the apparent strength obtained through Eq. (fJJ becomes 
exceedingly large. 

In conclusion, we have shown that while the phonon- 
induced electronic self-energy of graphene is isotropic and 
consequently the angular dependence of the electron- 
phonon coupling strength is negligible, the apparent 
electron-phonon coupling strength extracted from the ex- 
perimental angle-resolved photoemission spectra using 
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Wavevector (arbitrary unit) 

FIG. 4: (color online) Quasiparticle bandstructures of model 
systems including the electron-phonon interaction (dashed red 
lines). In (a) and (b) the bare electronic bands (dash-dotted 
blue lines) are assumed to be linear and quadratic, respec- 
tively. In each case, the actual electron-phonon coupling 
strength is set to A = 0.2. The horizontal solid lines represent 
the phonon energy u p h = 0.2 eV, and the energy AE — 0.3 eV 
below the Fermi level (Ef — 0) at which the slope is taken 
to calculate the apparent strength A app . The solid green line 
segments are tangential to the quasiparticle bandstructure at 
E = or E = —AE. 



Eq. ([4} exhibits significant anisotropy due to the curva- 
ture of the underlying bare electronic bands. Our analy- 
sis indicates that the band curvature is a crucial ingredi- 
ent for the interpretation of angle-resolved photoemission 
spectra. For example, the present result may carry impli- 
cations in the interpretation of many-body renormaliza- 
tion effects along the antinodal cuts in the photoemission 
spectra of cuprate superconductors, due to the presence 
of a saddle-point van Hove singularity along the Cu-0 
bond directions [34], HH • 
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